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Abstract 

Reactor neutrinos play an important role in determining parameter 613 in the lepton mixing 
(PMNS) matrix. A next important step on measuring PMNS matrix could be to build another 
reactor neutrino experiment, for example, DaYa Bay in China, to search the possible oscillations 
via sin^ 26'i3 and A77i^3 . We consider 4 different schemes for positions of three 8-ton detectors 
of this experiment, and simulate the results with respect to an array of assumed "true" values 
of physics parameters. Using three kinds of analysis methods, we suggest a best scheme for this 
experiment which is to place a detector 2200to ^ 2500m symmetrically away from two reactors, 
and to put the other two detectors closer to their corresponding reactors respectively, almost 
at a 100m ~ 200m distance. Moreover, with conservative assumption on the experimental 
technique, we construct series of allowed regions from our simulation results, and give detailed 
explanations therein. 
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I Introduction 



An next step in the exciting field of neutrino physics would be to improve current measurements 
and to measure some of the remaining unknown parameters in the full 3x3 leptonic flavor mixing, 
which is called the PMNS mixing matrix. There are important differences between the PMNS 
and quark CKM matrices. Which may be essential for our understanding the underlying physics. 
In addition to three masses mj, there are 6 free parameters in the matrix. We may parameterize 
UpMNS [1] as follows: 

UpMNS = Ul2 X U23 X Ui3 X Umajor- (1-1) 

COS 6*12 sin 012 
Ui2 = I -sin 6*12 cos 6*12 | (1.2) 
1 

1 

U23 = I cos 023 sin 023 I (1.3) 

- sin 023 cos 023 

Ui3 = \ 1 I (1.4) 

cos 013 




Umajor = \ 6"*" | (1.5) 

e-'^ 



U23 and the value of Am23 have been measured by the Super-Kamiokande [2] and K2K long 
base line experiments [3]; while U12 and the value of Amf2 are parameters of the confirmed solar 
neutrino MSW [4] solution. Umajor is the possible majorana phase matrix, its values and the overall 
mass scale will require kinematical and neutrino-less double-/? decay measurements. U13 is the next 



goal for the experiments since it will tell us, in particular, if CP violation is possible in the lepton 
sector. We can expect it from the very long base line experiments such as JPARC, H2B [5, 6] and, 
new generation of reactor neutrino experiments. 

The CHOOZ experiment [7, 8] only gives an upper limit on the mixing angle ^13 (sin^ 2^13 < 
0.10 at 90% confidence level). There are attempts to find this mixing angle in LBL accelerator 
experiments [5, 6], or in three neutrino analysis of solar neutrino [6, 9, 10], but the precision is very 
difficult to achieve. Since Am^2 ^ ^"^Is' must happen that Am^g ~ Amgg. A new generation 
of reactor experiments has been proposed to search for i>e disappearance at baselines of 1 2 km 
corresponding to this value of Am^. To improve on the mixing angle sensitivity achieved by Palo 
Verde and CHOOZ, proposals for reactor ^13 experiments include a large detector to reduce the 
statistical error, and also a second detector positioned very close (~ 100 m) to the reactor. The 
near detector would precisely measure the incident flux, providing to drop out many systematic 
uncertainties in the flux calculation. This also requires the detectors to be made identical and/or 
movable. Sensitivity down to sin^ 26*13 ~ 10^^ seems within grasp. Such experiments were discussed 
in literature [11, 12]. A practical possibility is a reactor experiment at DaYa-Bay, which is located 
near a special economic zone in Guang-Dong Province in southern China. There are nuclear power 
plants in that area. 

The knowledge we have about neutrino mixing is powerful to judge Grand Unifled models, such 
as the most inspiring SO(IO) GUT. Starting from the lepton quark symmetry in this model, one 
is able to obtain a bi-large neutrino mixing pattern via see-saw mechanism. It leads to a non-zero 
sin^20i3, e.g., at about 0.09 in paper [13], which is out of CHOOZ's limit but is very easy to 
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discover in DaYa-Bay experiment within one year operation. 

In this paper we will describe the importance of reactor neutrinos in determination of sin^ 2^13 
- Am^g; we will concentrate on possible new China experiment and its goal - finding sin^ 2^13. The 
paper is organized as follows: first, we will describe a possible reactor experiment in DaYa-Bay (see 
fig. 1), which is a kind of upgraded CHOOZ experiment [11, 14, 15]; its possible systematic and 
statistical uncertainties are analogies to the CHOOZ's one. Next we will explain different methods 
of analysis, and importantly, our arrangement of the detectors' positions, accompanied by our 
simplified Monte Carlo simulations on the possible sensitivity regions for sin^2^i3; and discovery 
potentialities are discussed therein. Finally, we will give our conclusion of the paper. 

II Reactor Neutrinos 



In a reactor, anti-neutrinos are released by radioactive isotope fission; the total neutrino spectrum 
is a rather well understood function of the thermal power W, the amount of thermal power wt 
emitted during the fission of a given nucleus, and the isotopic composition of the reactor fuel fi, 

The index i of fi stands for 4 isotopes such as ^^^C/, ^^^[/, ^^^Pn and '^'^^Pu. The (dN/dE) is the 
energy spectrum of the fissionable isotope, it can be parameterized by the following expression [16] 
when > 2MeV: 
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the coefficients aj depend on the nature of the fissionable isotope. KaniLAND is a scintillator detec- 
tor, where electronic anti-neutrinos are detected by free protons via inverse /3— decay reaction[16], 

Ve+p^e'^+n (2.3) 

In the limit of infinite nucleon mass, the cross section of this reaction is given by cr{Eu) = kE^+P^+ , 

where E^+, are the positron energy and momentum respectively and k can be taken as A; = 

9.55 X 10~^^ cm? /MeV'^. The anti-neutrino events are characterized by the positron annihilation 

signal and the delayed neutron capture sign [17]. 

Prom the reactor to the detector, massive neutrinos oscillate on the way and change their flavor 

composition to a certain extent. The anti-neutrino Vf. can oscillate to other flavors via Amfg and 

Am^g. In principle, reactor neutrinos are able to give us information about Am^g & Am^g and 

two mixing angles sin^ 20i3 &: sin^20i2, which are almost all the neutrino oscillation parameters 

except sin^ 2^23 • This advantage of reactor neutrinos is due to its low energy, which is several MeV, 

comparing with accelerator neutrinos' GeVs. Reactor experiment like DaYa Bay is sensitive to 

sin^ 2^13 and Amfg. This is the second order oscillating effect of the reactor neutrinos, comparing 

with the oscillations induced by 612 and Amfg- Since the energy of detectable reactor neutrino is 

about MeVs, from atmospheric oscillation experiments, Amfg is in the range of 1.3 ~ 5.0 x 10^'^ eV^ 

at 99%C.L. [2], with the best fitted point 2.5 x 10~^eV^, so the expected maximum oscillation is 

reached at about 1500 m. Then the survival probability is reduced to the two flavor neutrino case: 

, 9 N 9, . s 9 /1.27Am?o(eV2)L(in)\ , , 

L, ^13, Amfg) = 1 - sin2(20i3) sin^ _ilL_LAi j (2.4) 

where sin^ 2^13 is the only unknown parameter here, which is our major interest in this paper. 



5 



Ill Possible Reactor Neutrino Experiment in China 



DaYa-Bay is one of the Chinese running reactor power plants, located in Guang-Dong province, 
south China, and it is quite close to Hong Kong. The DaYa Bay nuclear power plant consists of 
two twin reactor cores, separated about 1200 m from each other, one is called DaYa, the other is 
called LingAo, each core can generate a thermal power of 2.9 GW, making a total of 11.6 GW, and 
a third twin core is planned to be on line in 2010. The reactors are built near a mountain, so it is 
possible to build a near experiment hall with an overburden of 400 mwe at a distance of about 300 
m to the core and a far hall with an overburden of 1200 mwe at a distance of about 1500 ~ 2000 
m to the core. This is a certain improvement in comparison with the CHOOZ experiment which 
has a value of 8.5 GW in thermal power, and 300 mwe in rock overburden. 

Three liquid scintillation calorimeter detectors can be placed from distances between several 
hundreds meters (near detectors) to about 2 km (far detectors) from the cores as shown in figs. 4 
- 7. High intensity and purity electron anti-neutrino flux from the reactor core is detected via the 
inverse beta-decay reaction eq. (2.3), the signature is a delayed coincidence between the prompt 
e+ signal and the signal from the neutron capture, the target material can be the Hydrogen-rich 
(free protons) paraffin-based liquid scintillator loaded with Gadolinium, which is chosen due to its 
large neutron capture cross section and to the high 7-ray energy released after n— capture. The 
rock overburdens are able to reduce the external cosmic ray muon flux by a fact of more than 
300. To protect the detector from natural radioactivity of the rock, the steel vessel should be 
surrounded by ~ 100cm of low radioactivity sand and covered by ~ 15cm of cast iron. There are 
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three concentric regions in a detector: as shown in fig. 2: one is a central 8 tons working target 
in a transparent Plexiglas container fihed with a 0.09% Gd-loaded scintillator; the second should 
be an intermediate region, equipped with hundreds eight-inch PMT's, used to protect the target 
from PMT radioactivity and to contain the gamma rays from neutron capture; the third can be an 
outer optically separated active cosmic-ray muon veto shield equipped with two rings of tens eight- 
inch PMT's. All these arrangements could significantly decrease the most dangerous background. 
Moreover, the particular coincidence between the prompt positron signal and the delayed neutron 
capture signal could significantly reduce the accidental background. Rates from background can 
be suppressed down to 0.2 ~ 0.3 events per day per ton, which corresponds to more than 60% of a 
detection efficiency. 

Besides the uncertainty in the efficiency, there are other uncertainties contained in the ex- 
periment: the reaction cross section uncertainty from an overall and conservative uncertainty on 
integral neutrino rate, 1.9%; the number of target protons uncertainty (0.8%), mainly because of 
the difficulty in determining the hydrogen content; the overall precision on the thermal power is 
claimed to be 0.7%; the uncertainty of the average energy released per fission of the main fissile 
isotopes is 0.6%. The five uncertainties presented above are all overall effective, and we could com- 
bine them to a total effect with a value of 2.7%. However, as the uncertainty of the energy scale 
(1.1%) affects the experimental result varying with respect to the energy bins, we have to specially 
consider its infiuence. 

It is supposed that only one of the two nucleon power plants would run for the most of the 
total time of the experiment. This allow us to know from which reactor the signal is generated. 
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Otherwise, the detector should be desired to be able to distinguish the events induced by the first 
reactor binary-core from the other. In the paper [18], a good determination of the anti-neutrino 
incoming direction is discussed. It is based on the neutron boost in the forward direction via the 
inverse beta-decay reaction, which is induced by the incident neutrino; the kinetic energy of the 
neutron remains even after collisions with protons inside the detector. 

IV Three methods of Analysis 

Let us suppose that the experiment like DaYa-Bay will measure the positron spectrum in 7 bins 
(from El to Ej). For a mean reactor-detector distance L^, the rates can be written as 



Sk{E,Lk,9,Am^) = J h{L,Lk) j a{E^)S{E^)P{E^,L,9,Am^)r{E^+,E)e{E,+ )dE^dE^+dL 

(4.1) 



The fissile isotope composition varies with respect to the working time, because the 4 isotopes 
burn a bit differently. The expected number of neutrino is a function of working time, reactor 



where 
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power and a constant background. 

Ni = (B^ + W,, ■ y/, + W2^ ■ Yl) . AU (4.2) 

where the index i labels the run number or the working time information; Atj is the corresponding 
live time interval; j labels the detector number; is the background rate, which is assumed to 
be a constant with time in our study; {Wu, W2i) are the thermal powers of the two binary-core 
reactors in GW and iYu, ¥2%) are the positron yields per GW per day induced by each reactor. 
Since we are considering several detectors, it is convenient to factorize Y^^ and B^ into functions 
which separate the factors which are independent of the detectors' size: 

yL = {^ + rf,d-^^^-Xi-r\ {A.2.) 

Bi = y>-T^, (4.4) 

where A; = 1,2 is the index of the reactors, associated with tons of available detector material, 
?7^j stands for corrections from the reactor's differences on the fissile isotope composition and the 
positron efficiency correction. The complications varying with time are all involved in 77-^^; this 
effect will be considered in real measurement and will not be taken into account in our simulations. 

represent the positron events contribution in 1 ton fresh detective material from a fresh reactor 
located 1km away, with its thermal power equal to IGW within one day. In the following we 
compare the two positron spectra X-j^ which are the measured one and the oscillation expected one. 
It can be parameterized by separating the oscillation term from the no-oscillation one: 

xiiEi,Li,e,Amj,) = Xo{Ei)P{Ei,Li,e,Aml,), (4.5) 
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where Xq{Ei) is the no-oscillation positron spectrum, unitary for all the detectors in our definition, 
and El label the visible energy bins: 

Xo{Ei) = ^^^^l^^2 I <y{E,)S{E,) I r{E,+ ,E)e{E,+ )dE,+ dE, (4.6) 

and P is the survival probability averaged over the energy bins and the finite sizes of both the 
detector and the reactor core, 

i . 2. [ Sk(E,Li,6,Am^)dE , ^ 

P(Ei,Li,0,Am^) = ^ ^ ' ^ 4.7 

^ ' ^ fSk{E,LlO,0)dE ^ ' 

With these definitions, we can begin to investigate our detector systems' power. In the reactor 
neutrino oscillation experiment like DaYa-Bay, we suggest that several detectors located at different 
places have no difference between them, which means they have a same design; thus the systematic 
uncertainties are supposed to be the same in order to cancel out some negative effects by comparison. 
In this paper, without losing generality, we suppose that: three detectors have 8 tons of Gd-loaded 
scintillator each; two nuclear reactors each working at 5.8GW of its full thermal power. 

To test an oscillation hypothesis (Am^, s\v?26) in our experiment, we construct a function 
including the 6 positron spectra measured and the oscillation expected in a 42-element X array, as 
follows: 

X = {Xl,X'l,Xl,X\,Xl,Xl), (4.8) 

x| = (X/(i^0,^(i?2),X/(S3),X/(i54),X^(i?5),^/(i^6),^/(^7)), (4.9) 

where the subscript i is the reactor's number, and the superscript j is the detector's number. Com- 
bining the statistical variances with the systematic uncertainties related to the neutrino spectra. 
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the 42 X 42 covariance matrix can be written in a compact form as follows: 

Vij = 5^^j{cj} + + ((5,,,_2i + 5ij+2i) • trS, J = 1, • ' -42) (4.10) 

where Uj are the statistical errors, Ui are the corresponding systematic uncertainties, and <t^2 
the statistical covariance of the reactor 1 and 2 yield contributions to the i-th measurement. 

A. 

There are two kinds of uncertainties influencing all the measurements which can't be ignored, the 
overall normalization uncertainty which is a a ~ 2.7%, and the spectrum affective energy-scale 
calibration uncertainty which is ag = 1.1%. We consider these two as in eq. (4.11) to reach the 
minimum for a set of proposed oscillation parameters: 
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^(X' -a-X {gE, 6, Am^)) (yi)^^^. {X^ -a-X {gE, 0, Am^)) 

(4.11) 



where X^ stands for measured or simulated one, and X {gE, 9, Am^) for oscillation expected one. 
Here we have taken into account the statistic uncertainty, so we add up a hat on X*. The statistic 
error is the square root of the event number which depends on the detector size, the reactor power 
and the experiment's life time. This method uses all the experimental information available, and 
directly depends on the correct determination of the systematic uncertainties. Such uncertainties 
could be reduced by measuring the positron spectrum with a detector near the reactors. When a 
product of the reactor power, the detector size and the working time is large enough to reduce the 
statistical uncertainty, we could constrain this overall normalization coefficient and the energy-scale 
calibration precisely. 
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B. 

We compare the ratio of the two reactors' contribution in another way. Since the expected spectra 
are the same for both reactors in the case of no-oscihation, the ratio reduces to the ratio of the 
average survival probabihties. Detectors are assumed to be the same, the systematic uncertainties 
will be canceled out or almost canceled out and the remaining uncertainty we should take into 
account is just the statistical uncertainty. So we can construct x'b 

is a good quantity to judge an oscillation result from no-oscillation one in the case where the 
distances to every reactor for a detector are quite different. In the CHOOZ experiment the distance 
difference is just 116.7 m, so that the Xb ^^^^ powerful than Xa [''']■ '-'^'^ experiment design, 
the near detectors are expected to be placed at positions which give a larger distance difference; 
thus we can expect x\ is more powerful. 

c. 

There is an intermediate analysis approach between the and x^- It uses the shape of the positron 
spectrum, while leaving the absolute normalization free. Similar to approach A, this approach fits 
the two uncertainties, but leaving the normalization parameter unconstrained, that is: 



'XniOi Am^) = min 



^{X^ - a ■ TigE, e, Am^)) {V"')^^^^ {P -a-X {gE, 9, Am')) + (^^) ' 
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In the following, we will use these three analysis methods to constrain the oscillation parameters 
for the possible reactors and detectors systems. With these results, we are able to judge which kind 
of system and the corresponding approach is the most powerful one. 

V Examining Arrangements and Simulation Results 

Comparing with the CHOOZ experiment, DaYa-Bay has more powerful reactors, longer working 
time, more larger detectors. Let us suppose the nuclear power plants work at its full thermal power 
5.8(GW), the detector's available mass is 8 tons, and the experiment life time is set to be one year 
in our first simulation. We assume we had known the signal's corresponding reactor individually. 

1) In our first scheme, three detectors with distances from the reactors in the range of 400m ~ 
1700m are arranged asymmetrically as plotted in fig. 4. Every detector can give two different 
distances with its corresponding oscillation spectra. This scheme is based on an idea: the allowed 
region of Am^ (= Am^^^) given by the atmospheric data is still too big to determine the best 
position for all detectors; in order to consider all the possible Am^, we try to place the detectors 
at different positions which can provide more numbers of different distances. We test our system 
with some possible parameters chosen from the CHOOZ 's indistinguishable region; the points in 
that parameter region are shown in fig. 3. We simulate the experiment by using the values of those 
points one by one as the values of true physics parameters. Like in previous discussions, we assume 
each possible oscillation hypothesis to construct the excluded/allowed region at different confidence 
levels [19]. 

The Monte Carlo simulation results are presented in figs. 8, 9, using three analysis approaches 
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for the six representational possible "true value" points as signed in fig. 3. For our first detectors' 
arrangement, we can see that approach B is the most powerful one while A is less sensitive and 
C is the lowest. In the following discussion, we will not present our approach C results since A is 
similar and better. Using analysis approach B in First scheme, the experiment like DaYa-Bay is 
able to find an oscillation result for sin 2^13 bigger than 0.05 (at more than 90% C.L. ), if Am^g 
is at the region of ~ 2.5 x lO^^eV^; and the oscillation parameters can be even constrained to a 
small allowed region (fig. 9) for bigger sin^ 2^13 (as the last three points we simulated with). It is 
interesting that one can bound sin^ 2^13 to be bigger than 0.05 at 95% C.L., if the true parameter is 
located near the point (sin^ 2^13 = 0.1, Am^ = 2.0 x 10~'^eV^); and the best restriction on sin^ 2^13 
is reached when the simulated parameters is located near (sin^ 2(9i3 = 0.1, Am^ = 2.5 x lO-^eV^); 
there we could use this system combined with analysis B to constrain the oscillation mixing angle 
and the mass square difference to a precise region as we present in the subplot (3rd row, 2nd 
column) named as 'FirstBS' of fig. 9. 

2) In the second scheme[20], fig. 5, the three detectors are arranged symmetrically from the 
reactors, one is in the middle of two reactors, the other two are superposed on the perpendicular 
bisector equidistant 1500m from the two reactors. Since the two reactors are symmetric to every 
detector, approach B is disabled. Approach A can also distinguish sin^ 2^13 bigger than 0.05 if 
Am23 is near 2.5 x 10~^eV'^; but it can only constrain sm^20i3 at 95% C.L. and the maximum 
mixing is not excluded at 99% C.L. even at the most sensitive oscillation points, if by the DaYa-Bay 
experiment itself. 

3) In the third scheme [20], fig. 6, detectors are also disposed symmetrically, but from the near 
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detectors, the two reactors are not equidistant, so the analysis approach B is partly available. For 
one year's operation of the experiment, this scheme shows a very clear oscillation signal for sin^ 2^13 
larger than 0.05 at 99.73% C.L. (2nd row, 3rd column in fig. 10). 

4) Scheme three is quite good since it gives a higher confidence level than First and Second 
schemes if there is a positive result for a non zero sin^ 2^13. However it is still possible to improve 
it. As a conclusion, in this paper, for detectors' locations, we suggest a possible setting for the 
experiment like DaYa-Bay, which is called the Fourth scheme: an extension of the Third scheme, 
which is to put a 8-ton detector 2200m ~ 2500m symmetrically away from two reactors; and put 
the other two 8-ton detectors more close to their corresponding reactors respectively, almost at a 
100m ^ 200m distance; they are located on the line between the two reactors (see fig. 7). The 
reason for a 2200m ~ 2500m-detector is based on the most sensitive oscillation zone with respect 
to the range of present Am|]^, taking into account the whole energy spectrum effect of reactor 
neutrinos. The best way for the two near detectors to easily distinguish which reactor a neutrino 
signal comes from, is to put the other two detectors on the inner line between two reactors; thus 
two different neutrino sources are from two opposite directions. We use the Xb inethod to analyze 
Monte Carlo results, in which the far detector's data is not taken into account since this detector 
is symmetric to two reactors; this doesn't affect our major statement about the discovery potential 
of this experiment. However, the far detector is important in a real data analysis: using the single 
x\ method it will exclude big-sin^ 2^13 region, and give a precise allowed region in a combined 
analysis of Xa ^^'^ x'b'j it can also be used to implement much smaller systematic errors in other 
different analysis method not discussed in this paper. During 3 years of running, the experiment 
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like DaYa-Bay is able to discover a non zero ^13 if sin^ 2^13 > 0.02 at 95% C.L. (fig. 11), and at 
99.73% C.L. for sin^ 2^13 > 0.03. This result is able to exclude the most general SO(IO) GUTs, an 
inspiring Grand Unification candidate, if nature doesn't choose it. 

As we have seen, the experiment like DaYa-Bay described in this paper is simple, has no 
technical difficulty, and could have been realized even several years ago; but with the optimization 
to the detectors position, we could get a more precise result (sin^ 2^13 ~ 0.02) than what we ever 
had [7]. Moreover, it is possible that some improvements in the technique, more detectors and 
advanced methods could be used [20] to have more precise results. It is promising that, to a greater 
extent, this experiment could reach a precision of sin^ 2^13 ~ 0.01. If the detectors are constructed 
as movable objects, they can be used to measure the solar neutrino parameters (sin^ 2^12, Amf2) 
in the next phase. 

VI Conclusion 

The reactor neutrino experiment like DaYa-Bay offers an opportunity to discover a non zero ^13, 
another crucial step in particle physics after solving the solar neutrino problem. We arrange four 
schemes for the three 8-ton detectors' locations, and select the fourth scheme as our suggestion 
for the experiment. In the First scheme, with respect to two reactors, we place three detectors as 
asymmetric as possible in the distance range of 400m to 1700m, in order to have more oscillation 
distances. We relax a systematic uncertainty to totally a few percent (much bigger than one percent 
that is a possible but difficult achievement), which is already reached by CHOOZ's technology. 
During three years of data taking, the simulation result shows that a discovery ability of this 
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scheme is sin^ 2^13 > 0.03; while one year operation can give for sin^ 2^13 a hmit of 0.05. The 
Second scheme is to place two 8-ton detectors at ISOOm in the same place, while a third 500m 
symmetrically away from the two reactors. The Third scheme is to put a 8-ton detector 1500m 
symmetrically away from the reactors; for the other two, both of them are 300m away from one 
reactor and 1237m from the other reactor, as shown in fig. 6. Both of these schemes are able to 
reach a limit of sin^ 2^13 at 0.05, during one year of data taking; moreover, the Third scheme gets 
this sensitivity with the highest confidence level of 3cr (fig. 10, 2nd row, 3rd column). We conclude 
that for a discovery potential, the Third scheme is a bit better; for a precise measurement after 
discovering a non- vanishing ^13, the First scheme is better. Furthermore, we suggest as the best 
possible location of detectors for the experiment like DaYa-Bay, the Fourth scheme: an extension 
of the Third scheme, which is to put a 8-ton detector 2200m ~ 2500?n symmetrically away from 
two reactors; and put the other two 8-ton detectors more close to their corresponding reactors 
respectively, at about 100m ~ 200m distance; and they are located on the line between the two 
reactors. The Fourth scheme will be able to discover a sin^ 2^13 > 0.02 at 2o" level (fig. 11, 1st 
row, 2nd column), for 3 years of running the experiment like DaYa-Bay. With improvement in 
technology and better budget, the sensitivity to sin^ 2^13 can be even better. 
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VII Appendix: The minimization method in approaches Xa 

Xc 

In section IV, we defined in Eq. (4.11) that the Xa minimum of the fitting of a,g for the 

oscillation parameters. Because the energy-scale calibration factor g is involved in the integral for 
X, such as in Eq. (4.6), it is troublesome to get an analytical expression for the minimization, and 
the numerical computation is also insufferable if scanning the parameters plane is necessary. In 
order to get a precise minimum quickly, we assume that: 

XH,il.E,lA^^) ^ = /(«).(.-!) (7,1) 

where f{Ei) is the ratio of the two differences, gained from the no-oscillation case, since the oscil- 
lation effect is very small even though it can be detected in our system, these ratios are suitable for 
the oscillation case. We check the linear assumption at every energy bin, it holds when g changes 
near 1.0 in the range of 5%. This property is good enough when the uncertainty of g is just 1.1%. 
For convenience, we rewrite Eq. (4.11) as 

xl = inmWia,g) (7.2) 

where 



W{a,g) = J2 {(X' - aT{gE)) {V-')^^^ {X^ - aX\gE))) + (^X + (^X (7.3) 

Here we have omitted the oscillation parameters in the bracket, keeping in mind that X followed by 
bracket is the oscillation prediction, while single X is the experimental or Monte Carlo simulated 
result. Using the linear assumption, we can get W[a,gys partial derivatives daW{a^g),dgW{a,g) 
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easily, moreover, the numerical computation is just the multiplication of vectors and matrices. 



dW{a,g) 
da 



2 J] [-X\g = 1) {V~^)^^^X^ -{g- l){fXy {V-').^^X^ 



+aX\g = 1) {V-')^^^X^{g = 1) + 2a{g - l){fXy {V"') ^^^^ X^ {g = 1) 



+aig-if{fxy {v~')^^^{fxy 



a - 1 



(7.4) 



dW{a,g) 
dg 



-2 [a{fxy x^ - a\fxy x{g = ly 



-a\g - l){fXy (y-i) . ^. ifXy] - 2^ (7.5) 



'9> 

where {fXy is f{Ei)-X'^{g = l,Ei). Starting from (a = I, g = I), diivenhy {-daW{a, g), -dgW{a, g)) 
at this point, and using appropriate iterative step length, the minimization of W{a, g) for an oscil- 
lation hypothesis is quickly achieved. 
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Figure 1: The picture of the DaYa-Bay reactor. 
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Figure 2: Structures of a detector for reactor neutrino experiment, the experiment like DaYa-Bay 
can have three such 8-ton detectors. The picture is from CHOOZ detector: "region 1" contains 
5-ton target material, "region 2" protects the target from PMT radioactivity and contains the 7-ray 
from neutron capture, "region 3" is used to separate active cosmic-ray muon veto. This figure is 
taken from paper [7]. 



23 



analysis A 

analysis E 

— analysis C 



10 



0.1 0.2 



0.3 



0.+ 0.5 



0.6 0.7 0.S 0.9 



sin '(26) 



Figure 3: Six points whose coordinates are considered as real physics parameters for the first 3 
schemes' simulations, in CHOOZ's excluded region. First ponit represents no-oscillation, labelled 
as "0". The other five points below will cause oscillation effect; their parameters (sin^2^i3, Am? 
[eV2] ) are 1^(0.02, 2.0 x IQ-^), 2^(0.05, 2.5 x IQ-^), 3^(0.1, 2.0 x lO'^), 4^(0.2, 1.7 x lO'^), 
5^(0.1, 2.5 X 10-3). 
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Figure 4: First scheme: we arrange three detectors in distances from the reactors in the range of 
400?7i ~ 1700m unsymmetrically. 
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Figure 5: Second experiment scheme: three detectors are placed symmetrically according to two 
reactors, one is on the middle points, the other two are superposed on the perpendicular bisector 
equidistant 1500m from the two reactors. 
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Figure 6: Third scheme: we arrange three detectors symmetrically also, but from the near detectors, 
the two reactors are not equidistant. 
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Figure 7: Fourth scheme: similar to scheme 3, but two near detectors are put on the hne between 
two reactors; and the distance of the far detector is enlarged to 2200 ~ 2500 m. 
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Figure 8: The allowed regions at different confidence levels for the first scheme. The labels "First", 
"Second", "Third" stand for different schemes; "A", "B", "C" for different analysis methods, and 
"0 ~ 2" for the simulation input parameters with their value presented in the figure, in a plot, it 
is denoted by a plus sign. From this point to the exclusion area, the confidence levels of the four 
regions are 90%, 95%, 99%, 99.73%. 
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Figure 9: Continue plots of fig. 8 
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Figure 10: The allowed regions (shadows) of our simulation for second and third schemes, same 
notations as previous plots. 
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Figure 11: Similar to previous figures, these are results for the Fourth scheme with 3 years' exper- 
iment. The title "FourthBOx" stands for Xb analysis of the simulation with the parameters in the 
figure. 
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